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Abstract 

Using methods from effective field theory, we have recently developed a novel, systematic 
framework for the calculation of the cross sections for electroweak gauge-boson produc- 
tion at small and very small transverse momentum qt, in which large logarithms of the 
scale ratio my /(It are resummed to all orders. This formalism is applied to the produc- 
tion of Higgs bosons in gluon fusion at the LHC. The production cross section receives 
logarithmically enhanced corrections from two sources: the running of the hard matching 
coefficient and the collinear factorization anomaly. The anomaly leads to the dynamical 
generation of a non-perturbative scale ~ 777,^ g-const/os(mH) ~ 8GeV, which protects 
the process from receiving large long-distance hadronic contributions. We present de- 
tailed numerical predictions for the transverse-momentum spectrum of the Higgs boson, 
finding that it is quite insensitive to hadronic effects. 



1 Introduction 



The transverse-momentum spectra of electroweak bosons are among the most basic observables 
at hadron colhders. At large transverse momentum qt these spectra can be computed in fixed- 
order perturbation theory. On the other hand, if the transverse momentum is much smaller 
than the boson mass, higher-order corrections are enhanced by large Sudakov logarithms and 
must be resummed. This resummation was achieved a long time ago by Collins, Soper, and 
Sterman (CSS) [l| and has been implemented to high accuracy in several numerical codes 
jij-jij. For vector-boson production, we have revisited the resummation in the context of soft- 
collinear effective theory (SCET) and have derived an all-order factorization theorem for the 
cross section at small qx 0. The factorization theorem for the differential cross section is 
affected by a coUinear factorization anomaly, which generates an additional dependence on 
the hard momentum transfer in the low-energy theory. It was shown that in position space 
this extra dependence takes the form of a pure power of the gauge-boson mass. Relating our 
result to the traditional CSS formula then allowed us to determine the last missing ingredient 
needed for resummation to next-to-next-to-leading logarithmic (NNLL) accuracy. 

In the recent paper we have used our formalism to analyze the differential cross section 
da/dqT at very small transverse momentum, where it exhibits quite remarkable properties. 
In this region, the spectrum is genuinely non-perturbative but dominated by short- distance 
physics and therefore calculable. For a boson of mass my, long-distance effects are suppressed 
by a dynamically generated scale ~ my e-const/as(mv)^ which is close to 2 GeV for the case 
of Z-boson production. While the underlying mechanism was identified a long time time ago 
[ij, our framework has allowed us to systematically compute corrections also in this region. 
In j3] we have performed a detailed phenomenological study of Z-production at the Tevatron 
and the LHC and have investigated the numerical impact of long-distance effects. 

The region of low transverse momentum is also important for the study of the Higgs boson 
and its properties, since the background is reduced when additional radiation is vetoed. In 
practice, this is done by imposing a jet veto. Several recent papers have considered resumma- 
tion for the Higgs-boson cross section in the presence of such a cut js iS], and the resummation 
is now known to NNLL accuracy 10|, ll^l • In the present paper, we extend the formalism of 

to the transverse-momentum spectrum of the Higgs boson and perform the resummation 
of the spectrum to NNLL accuracy. Much of the extension is straightforward and some of 
the necessary ingredients were already given in jof, but there are a few interesting differences 
to the vector-boson case. First of all, while there is a single way to combine the spins of the 
incoming quarks to produce a vector boson, there are two ways to combine the gluon spins to 
produce a spin-zero state. As a consequence, the cross section is not just given by the product 
of two collinear functions as in the usual CSS formula, but a sum of two products of coUinear 
functions describing the two production mechanisms [iJ]. We show that the collinear anomaly 
is the same for both structures, and that the dependence on the large scale rriH therefore arises 
as an overall factor in position space. We then compute the collinear functions at one-loop 
order. The other feature, which distinguishes Higgs production from the vector-boson case, is 
that the infrared protection mechanism discussed above is much more efficient. The numer- 
ical value of ~ 8 GeV is significantly higher than in the Z-boson case, and we show that 
long-distance hadronic effects have almost no impact on the Higgs-boson spectrum. We have 
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implemented our resummed results for Drell-Yan, W, Z, and Higgs production in a public 
code CuTe 15| and give plienomenological predictions based on this program. 



2 Factorization and resummation 

We consider the cross section for the production of a Higgs boson with mass mn and transverse 
momentum qr = \ci±\ in gluon fusion at the LHC The derivation of the factorization formula 
for the cross section proceeds exactly as in the case of the Higgs-pro duct ion cross section 
defined with a jet veto, which we have recently considered in [10|. Our analysis there has 
been performed at fixed q± and rapidity y of the Higgs boson, and the integration over the 
boson phase-space was carried out at the end. We can thus immediately use the result for the 



factorized cross section obtained in 10 , which reads 



da = ao{fi)C^{mlfi)\Cs{-ml,fif dy ^ I d V e"^^-^^ 



TS (27r)2 J ^-^^ 
X 28^" {^i, x^, /i) i3g^^(^2, x±, /i) S{x^, fi) , 
where ^1^2 = a/^c^^ and r = {mjj + |g^|)/s. The Born-level cross section is 

where y/s denotes the center-of-mass energy of the LHC and v is the Higgs vacuum expec- 
tation value. The Wilson coefficient Ct multiplies the effective ggH operator obtained after 
integrating out the heavy top quark, while the hard matching coefficient Cs arises when this 
operator is matched onto an effective two-gluon operator in SCET. Moreover, we have defined 

^r(e,a:x,/i) = / rfte-^«*^-^ {Pm^-'clitn + x^) |X,) {X^\X^m \P{p)) , 

= 5^ (0 I {SiSnY\x^) |X.) (X,| (4^n)'"(0) |0) . (3) 

Xs 

Here Ad. is the gauge-invariant effective gluon field of SCET, and Sn, Sn denote soft Wilson 
lines. The function B'^^ describes the structure of the jet of collinear particles inside one of the 
colliding protons (the one moving along the light-like direction n^), which is probed at small 
transverse distance x^. The corresponding function B'^^ for the second beam jet, consisting of 
anti-collinear particles (moving along n), is given by the same formula with the replacements 
rz — )■ n and c — ?■ c. The soft function S describes the physics of soft gluons emitted from the 



colliding beam particles. In [10[, the sum over hadronic intermediate states was restricted by 



the jet veto, while in the present case the sums in ([3]) are completely inclusive. 
2.1 Collinear anomaly 

The beam-jet functions iS^'^, B'^^ and the soft function S suffer from light-cone divergences, 
which are not regularized by the conventional dimensional regularization procedure. These 
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singularities cancel in the product of the three functions, but in order to make the individual 
objects well-defined an additional regularization is required j^. At the end of the calculation 
the regulator can be removed, but a non-trivial effect remains: The additional regularization 
breaks a rescaling symmetry exhibited at the classical level by the effective Lagrangian of 
SCET (i.e., the property that the coUinear matrix elements are invariant under a rescaling 
of the anti-collinear momenta and vice versa), which is not restored in the limit where the 
regulator is removed because it is spoiled by quantum effects. This "coUinear factorization 
anomaly" manifests itself through a dependence of the product of soft and beam functions on 
the Higgs mass - the large scale in the problem. For the Drell-Yan process, the all-order form 
of this anomaly was derived in j6|, and it was shown that the dependence on the hard scale 
(the electroweak gauge-boson mass in this case) takes the form of a pure power in x_[_ space. 
We will now extend the derivation to the case of Higgs production and show that the form of 
the anomaly as well as the corresponding anomalous dimensions are spin independent. 

The simplest way to introduce the additional regularization in such a way that gauge 
invariance and the factorization properties of the cross section are maintained is to regularize 



the phase-space integrals analytically by replacing [16 

j d'^k 5(P) ^(A;°) ^ J d'^k 5(P) ^(A;0) (^^^ (4) 

in the sum over intermediate states in ([3]). The scale v is inserted to restore the canonical 
dimension of the integrals, in analogy to the scale of dimensional regularization. Instead of 
using a single light-cone component k^ = k ■ n m the regulator term, one could also use the 
sum k^ + k_ = 2kQ. This last form is similar to what is used in the "rapidity regularization 
scheme" proposed in [17]. For pertubative computations, using a single light-cone component 
is simplest, since SCET Feynman diagrams typically already contain light-cone denominators 
involving k+ or A;„. With the form (jlj) of the regularization, one finds that the soft function 
is given by scaleless integrals, and thus S{x±,jj) = 1 to all orders in perturbation theory 0]. 
We will no longer write it explicitly in the rest of the paper. If one were to use the energy 
ko instead of /c+ in (jlj), the soft function would be non-trivial, but can be absorbed into the 
beam-jet functions without loss of generality. 

The light-cone component fc+ of the anti-collinear particles moving along the direction is 
large, ~ mn- Expanding in the regulator a, the dependence of the anti-collinear beam-jet 
function on the regulator scale u thus takes the form ln(i//miy). On the other hand, the fc+ 
component of the coUinear partons is small, fc+ ~ Qt/^h ~ ^/{xrTnH)-, where = —qj_ and 
\. In the collinear beam-jet function the dependence on the scale u thus arises in the 
form \n{i> x'^uih)- The requirement that the physical cross section ([1]) must be independent of 
the analytic regulator scale z/ can then be expressed as 

^^^r(^i>a;±,ln(z/x|mH),/i) i3r(6,a;±,ln-^,/ij =0. (5) 

In the factorization theorem ([1]) the Lorentz indices of the beam functions are contracted. 
The fact that u independence also holds without contracting the indices follows by consider- 
ing the factorization theorem for the production of a general color- neutral tensor field h^"^. 
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The corresponding factorization theorem has the same structure as ([T]), except that the hard 
matching coefficient \Cs\'^ would now depend on the Lorentz indices of the tensor fields in the 
initial and final states. Since the logarithms in (|5]) have different arguments, the cancellation 
of the V dependence among the different factors imposes a non-trivial constraint on the mu 
dependence of the product. As explained in detail in 0, 18|, the above equation implies that 
the dependence of the product of the two functions on uih must be power like. We can thus 
rewrite the product in the form 



H 



ft2 



rriH 



(6) 



with Bq = 2e ""^ . The new beam-jet function /i) and the anomaly exponent 



, /i) are independent of rriH- 



Having determined the form of the anomaly, we now derive the scale dependences of the 
function /i) and the exponent Fgg{x'^,fi). Their anomalous dimensions can be in- 

ferred from the requirement that the cross section must be independent of the renormalization 
scale /i, which implies that the fi dependence of the product of beam functions must cancel 
against that of the hard function o'oCf\Cs\'^ in ([T]). This leads to the renormalization-group 
(RG) equations 



d 



dliafi 



dFgg{xl,H) 

d\n^ 

B^^{i,x^,li) 



2T 



cusp 



cusp 



[as) In 



27S 



(7) 



where T^gp is the cusp anomalous dimension in the adjoint representation, and 7^ is the 



anomalous dimension of the collinear gluon field as defined in [19l|. The fact that each com- 
ponent of Bj^'^{C,, x±, /i) renormalizes in the same way follows after considering the production 
of a general tensor field h'^'^ and using that the anomalous dimension of the hard function is 
spin- independent ji^]. Three-loop expressions for Fj^^p and 7^ can be found in 19 



Lorentz invariance implies that the renormalized beam-jet functions can be decomposed 



as 



9T 



1/' /-yi/^ /yj I' 



V2\ 



9T 



X , 



where the coefficients only depend on the invariant = —x^. It follows that 
25f(ei,xx,/i)5,^.(6,a;±,/i)= Yl Bf\i,,xl,ii)B^^\i,,xl,ii) 



(8) 



(9) 



n=l,2 



The presence of the two different tensor structures ([8]) contributing to the Higgs cross section 



was ffist pointed out in [IJ]. As a result, the Higgs-production cross section does not simply 



factor into a product of two beam-jet functions, but instead it involves the above sum of 



4 



products, which arises because the spin-0 state of two gluons is an entangled state. Rewriting 
the cross section ([1]) in terms of the new beam-jet functions, we obtain 



dq^ dy 



An 



-Fgg{x^,ll) 



(10) 



X 



which is analogous to eq. (18) of our paper [6| on the Drell-Yan cross section for the production 
of electroweak gauge bosons. In the above formula, the three disparate scales rrit 3> niH 3> x^^ 
appear in factorized form, and it will be possible to resum logarithms of their ratios by 
controlling the /x dependence of the various functions in (ITU]) using RG equations. After 
performing the Fourier integral the scales mu and qx get intertwined in a complicated way. 
At small qT this gives rise to interesting phenomena, such as a non-perturbative, but short- 
distance dominated dependence of the cross section on qr/mH and Osiqr) [gI. ItI. EiI| . 



2.2 Refactorization 

As long as the transverse displacement xt is much smaller than the scale of long-distance 
interactions in QCD {xt -C Aq^j-,), formula ffTOl) for the cross section can be simplified further, 

by matching the beam functions Sg"^ onto ordinary parton distribution functions (PDFs), 
thereby computing their dependence on in perturbation theory. The relevant matching 
relations read 

5f(e,4,/^)= E /'y/i2.(^,4,/^)0./p(e/^,/^), (n) 

i=9,<],9 ^ 

which is valid up to hadronic corrections suppressed by powers of Aq^.^, x'^. The cross section 
can then be written in the final form 



^=ao(/.)c,v?,/^)M-<,/.)r E L^j/-^ 

X Cgg^ij{zi, Z2,qT,mH, n) (t)i/p{ii/ zi,fi) 0j7p(6/ Z2, /i) 



;i2) 



where 



2 ^2 \ -F99iL±,as) 



(13) 



X 

n=l,2 



With a slight abuse of notation, we have traded the variables x^ and /i in the functions 



and for new variables 



in) 

■ iui licvv vaiiduicD 

I 



99 



...^ (14) 



without changing the names of these functions. This notation will be convenient for our dis- 
cussion below, and it conforms with the notation used in . Integrating the double differential 
cross section (|T2|) over rapidity, we find 



da 

dq^ 



ao{^i)Cj{ml,^)\Cs{-m]j,^)\^ I — Cgg^ij{z,ql,m%, ^i) ffij{T / z, ^i) , (15) 



where the parton luminosities and new kernel functions are defined as 

du 



u 



(j)^/N,{u,f^)4>j/N2{z/u,fl) 



^ du 



(16) 



^gg^iji'^iQ.Ti'^Hi lA — I Cgg^ij{u,z/u,qrp,mfj,n) . 



u 



The factorized cross sections (fT2l) and (fTSi) receive power corrections in the two small quantities 
q^/m'jj and AqQ^xf^, which will not be indicated explicitly in our equations. 

A dependence on the hard scale mn enters formula (|T2|) for the double-differential cross 
section in two places: via the hard matching coefficient Cs and via an x-r-dependent power of 
rriH under the Fourier integral in fll3p . The latter effect is due to the coUinear factorization 
anomaly As long as x'^ <^ ^qcd' anomalous exponent Fgg can be calculated in 
perturbation theory, and at least up to three-loop order it is related to the corresponding 
exponent Fgg appearing in the Drell-Yan case by the Casimir-scaling relation jij 



FggiLj 



Ca Cf 
Using the known expression for Fgg, we then find 

Fgg{L^,^^)=V^ 



(17) 



a,L^ + ai{ (3o^ + KL^ + da 1 + 



where Tq = AC a and /3o = y Ca — | Tpn j are the one- loop coefficients of the cusp anomalous 
dimension and /9 function, and 

[67 7r2\ 20 /202 \ 56 

contain the relevant two-loop information. Because of Casimir scaling, these coefficients take 
the same values as in the case of Drell-Yan production. 



At one-loop order, the kernel functions IgJ^^{z, L±, as) are given by 



ll^li{z,L±,as) = 5(1 - z)6gi 



1 + as 



9<- 



+ as 



r{2) 



I„J_Az,L^,as) 



as Ti'g^iiz) 



(20) 
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where 7^ = -/3o, 



= 8Ca 



l-z) 



1 + 



+ z{l-z) 



+ 2l3o5{l-z) 



(21) 



are the one-loop DGLAP sphtting functions, and the remainder functions lZg^i{z) and lZ'g^^{z) 
are given by 



TT 



n'g^giz) = -aV2Ca 



1 - z 



l-z 



(22) 



The expression for was calculated in [10|, while the result for is new 



2.3 Resummation of Sudakov logarithms 

The resummation of large logarithms in the cross section f[T^ is accomplished by evolving 
the hard matching coefficients Ct and C5 to a scale /x at which the kernel functions Cgg^^ij 
in (fT3|) can be calculated using a controlled perturbative expansion. The solutions of the 



10 . For our numerical work we use 



corresponding RG equations are discussed in detail in 
relations (A. 2) and (A. 4) from this paper. For the Drell-Yan case, the proper choice of the 
factorization scale /z has been discussed in [3]. Since the resulting expressions for Higgs 
production are completely analogous, we will not repeat details of the derivations here but 
rather summarize the main physical insights and quote the final expressions. The most naive 



choice would be to set ~ inside the Fourier integral in f[T3l) . in which case Lx_ would be 
a small logarithm for any choice of xt- There are several disadvantages to such a treatment. 
First, since xt is integrated over all possible values, there would be no clear meaning to the 
scale 11 in terms of a characteristic scale of the process. Second, setting the scale under the 
integral means that the integration unavoidably hits the Landau pole of the running coupling, 
giving rise to ambiguities in the numerical results. In the spirit of effective field theory, the 
scale /i should correspond to a physical scale in the underlying factorization theorem. We will 
choose it in such a way that on average the xr-dependent logarithm is small, and denote 
the corresponding value by /i ~ {xT^}). 

Naively, one would expect that the transverse momentum qt and average transverse sep- 
aration (xt) are conjugate variables satisfying qx ~ {'^t^)- While this is sometimes true, the 
general situation turns out to be more complicated. After integration over the factorized 
dependence on rriH and qx in f|T3l) gets intertwined in a complicated way, and this gives rise to 
the peculiar effect that the two scales qx and fi ~ (x^^) decouple for very small qx [7|. When 
this happens depends on the value of the coefficient 



In 



TT 



(23) 
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As long as < r/ < 1, one indeed finds that ~ (x^^) ~ qx, because contributions from large 
values xt Qt^ are suppressed due to the rapid oscillations of the phase factor of the Fourier 
integral, while contributions from small values xt <^ q^^ are phase-space suppressed. The 
situation changes, however, at very small transverse momentum, where with the prescription 
/X ~ gT the value of rj reaches 1. We denote by the value of fi where this happens, i.e. 



27r \ / 27r 



= niH exp ( — — ) ^ uih exp ( - , o \ ^ \ ) ' (24) 



where in the last step we have used the one-loop approximation for the running coupling. As 
long as is in the perturbative domain, one finds that at this scale (x^^) decouples from 
qx, and it remains a short- distance scale even in the extreme case where qx is taken to [3]. 
Changing variables from xt to Lj_ in the Fourier integral, one observes that the integrand 
exhibits a Gaussian peak with a width proportional to 1/ ^/a^. The condition that at the peak 
the logarithm L± = 0{1) implies that 1 — rj = 0{as), indicating that the factorization scale 
must be chosen in the vicinity of g*. We thus conclude that the proper scale choice is 

/i ~ {x^^) ~ max(gT, g^,) . (25) 

In our numerical work below, we will use /i = gr + g* as the default choice for the factorization 
scale. Solving the first equation in fl2^ numerically, we obtain q^ ^ 7.7GeV for mn = 
125 GeV, which is a short-distance scale well inside the perturbative domain. Due to the 
difference in color factors, this scale is significantly larger than in the case of the Drell-Yan 
production of electroweak gauge bosons, for which g* ^ 1.75 GeV [3]. 

It follows from these arguments that the transverse-momentum distribution of Higgs bosons 
is protected from long-distance physics even for arbitrarily small gy - a fact that in the 
context of the Drell-Yan process has been pointed out first a very long time ago in j^. The 
resummed perturbative series for the cross section generates the scale q^ dynamically, and even 
though this is a short-distance scale, it is related to the boson mass itih in a genuinely non- 
perturbative way. The scale g* also sets the magnitude of hadronic long-distance corrections, 
which turn out to be power-suppressed in the ratio Aqco/g*- We expect that these corrections 
are significantly smaller for Higgs production than for the Drell-Yan production of Z and W 
bosons. This expectation will be confirmed by our numerical studies presented below. 

The above discussion shows that we must distinguish two regions of transverse momenta. 
For qx ^ g*, the scale choice fi ^ qx prevents that the logarithms L± give rise to large 
perturbative corrections. It is then consistent to count these logarithms as L± ~ 1 and 
construct the perturbative series as a series in powers of a^. A different situation is encountered 
for gT ^ g*. Even though the scale choice /z ~ g* ensures that L± = 0{1) on average, the 
Gaussian weight factor allows for significant contributions to the Fourier integral over a range 
of larger values with a width proportional to 1/ ^/a^■ It is then necessary to reorganize 
the perturbative expansion by adopting the modified power counting L± ~ 1/ ^/a^. This 
implies that single-logarithmic terms {asL±)^ ~ a"^^ are always suppressed, whereas double- 
logarithmic terms {asL\)^ ~ 1 are unsuppressed and must be resummed to all orders. To 
keep track of this fact, we introduce an auxiliary expansion parameter e (which at the end is 
set to 1) and assign the power counting ~ e and L± ~ e"^/^. The terms contributing up 
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to 0{e) to the cross section have been derived in [7] using recursive solutions of the relevant 
RG equations. Adapting the resulting expression to the present case, we find that the hard- 
scattering kernels defined in (fT3|) can be written in the form 



Z2,qT,mH,fi) 



X 



dxT XT MxtQt) exp [gAijli L±, a^)] 



(26) 



n=l,2 



where 



as (2^1 + vK) + al [V^ + r//3o) /3o ^ 



(27) 



J el/2 



0(e3/2) . 



Note that we treat ln(m|^//i^) as a large logarithm and count rj defined in (12 3 p as an 0{1) 
variable. The auxiliary parameter e counts the order in resulting (for qx ^ q*) after the xt 
integral in (1261) has been performed. The two terms given in the first line are unsuppressed 
and must be kept in the exponent of the integrand in fl26|l . whereas the remaining terms can be 
expanded in powers of e^^^. It is important that the expansion is truncated at an integer power 
of e. The resulting integrals over the Bessel function can readily be evaluated numerically. 

We finally give the expressions for the coUinear kernel functions corresponding to our 
modified power counting. We find 



Ig^liiz, L±, as) = 6{1- z) 6gi - 



,1/2 



+ 



asTlg^iiz) + al (Vg^,{z) - 2(3oVl%{z] 



LI 



(28) 



0{e 



3/2> 



while J^^j coincides with /^^- in fl2Ul) up to higher-order terms in e. The corresponding 



contribution to 
quantities 



is of 0{e ) and can be neglected to the order we are working. The 



E 



3=9,1,1 



^ du 



u 



vfl,{u)vfUz/u) 



(29) 
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involve the convolutions of two DGLAP splitting functions. Following [^], we find 



^ 1 ' In^ ^-2(l + z)ln2 + 3(l-z) 



\-z 



z 



z 



+ 16C^/3o 



\-z 



z 1 — z 

+ + z(l - z) 



+ 32CFTFnf 



2(1 + z) \nz + l-z + 



3 z 

4 1-z^ 
3 z 



= 16CaCf 



1 + (1-^)2 



In 



2 + 6z + 3^2 



z 



+ 16C- 



1 + 



In 



'1-z 



z 

)2 21n^ 
- + + 2 



Inz - (1 
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si 
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3 



4^ 

T 



(30) 

Equations f l26|) - fl28|l are our main results. With the help of these expressions, large loga- 
rithms can be resummed at NNLL order for arbitrarily small transverse momenta. For larger 
qt values, the additional terms contained in (1271) and fl25]) reduce to higher-order terms pro- 



portional to and a^, which can be kept without doing any harm. Our formula thus provides 
a smooth interpolation between the regions of small and very small qt- In fact, it has been 
shown in [3| that the additional terms needed at very low qt serve an important purpose 
also at qr > q*, as they resum an asymptotically divergent, but Borel-summable subset of 
higher-order corrections in ag- 
in Figure dl we present predictions for the resummed transverse- momentum distributions of 
Higgs bosons with mass tjih = 125 GeV produced in gluon fusion at the LHC. The blue bands 
correspond to the results obtained at NLL order, in which case we keep terms up to 0{e^) in 
f l27|l and fl28|) and use leading-order (LO) expressions for the Wilson coefficients Ct and Cs in 
RG-improved perturbation theory. The green bands show results with NNLL accuracy, which 
are obtained by retaining all terms shown in (^7^ and f l25]) and using next-to-leading-order 
(NLO) expressions for Ct and Cs- We use v = 246.675 GeV for the Hig gs v acuum expectation 



23 . In our numerical 



value and work with MSTW2008NNLO PDFs with a,(Mz) = 0.1171 
results, we include finite top-quark mass effects by working with the exact leading-order cross 



section, which can be found e.g. in j2J]. Numerically, using rrit = 172.6 GeV, this leads to a 
6.6% increase in the cross section compared to the — )■ oo limit. We do not include the 
6-quark contribution and electroweak effects, which are of similar size but opposite sign and 
tend to cancel each other. The combined effect of these would be a reduction of the cross 
section by about 1.5%. To estimate the theoretical uncertainty of our predictions, we vary the 
factorization scale by a factor 2 about the default value fi = qr + q*- Our results also depend on 
the two hard matching scales fit and /i/j. At the scale /it the top quark is integrated out, giving 
rise to an effective ggH vertex. The associated matching corrections are small, and varying 
/jjt in the NNLL order cross section by a factor 2 about the default value fit = mt changes the 
result by less than 1%. The scale fih is associated with the hard momentum transfer in the 
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10 20 30 40 50 60 5 10 15 20 

qr [GeV] qr [GeV] 

Figure 1: Resummed predictions for the transverse- momentum distribution of Higgs bosons 
produced at the LHC at NLL (blue bands) and NNLL order (green bands). The factorization 
scale 11 is varied by a factor 2 about its default value ^ = qx + q*- The thick hues refer to the 
default scale choice. The discontinuities at qr = 1.8 GeV arise from the change from 5 to 4 
light flavors occurring when /i < ^{qx + q*) =4.75 GeV. 



ggH coupling, as described by the time-like gluon form factor. In |24j. |25|. we have shown that 
it is advantageous to evaluate the relevant matching coefficient Cs{—Tn'jj, Hh) with a time-like 
scale choice /i^ = —m'jj. This eliminates the large perturbative corrections arising when the 
gluon form factor is continued from space-like to time-like kinematics. When this is done, also 
the corrections to C si—m?^ , Hh) are of moderate size, and the effect of a variation of Hh by a 
factor 2 on the cross section is again below 1%. Since the variation of the factorization scale 
H leads to the largest scale uncertainties by far, we use it to generate the error bands in the 
plots, keeping the hard matching scales /it and fixed at their default values. We observe 
that, as expected, the results at NNLL order exhibit a significantly reduced scale dependence 
compared with those at NLL order. The NNLL corrections have a noticable effect and strongly 
enhance the cross section in the peak region. The left plot shows the spectrum da/dq^, while 
the plot on the right shows the distribution da/dq^, which tends to a constant at the origin. 
Note that at NNLL order the value of the intercept is predicted with very good accuracy. 
In 0], we have presented an explicit formula for the intercept for the case of the Drell-Yan 
process. 

We proceed to study the impact on long-distance hadronic effects on the transverse- 
momentum distribution, which for the case of Drell-Yan production of electroweak bosons 
are known to have a non-negligible impact 0,1,13. Following tS], we model these effects by 
noting that the beam-jet functions Bg^\^,x'^, n) in (fTO|) . which are nothing but transverse- 
position dependent PDFs, must vanish rapidly when the two gluon fields are separated by a 
transverse distance xt larger than the proton size. This motivates an ansatz of the form 

(e,x|,/i) = Aadr(xTANp)S(")P'='-*(e,4,/i), (31) 

where the perturbative functions Bg"'^^'^^^ carry all the scale dependence and are given by 
(ITTl) . whereas the hadronic form factor /hadr('") with /hadr(O) = 1 describes the fall-off at large 
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Figure 2: Comparison of the importance of long-distance hadronic effects on the differential 
cross sections da/dgr for Z-boson (left) and Higgs-boson production (right) at the LHC with 
y/s = 8 TeV. We adopt the Gaussian model (l32i) and vary Anp between and 1 GeV. The 
cross sections for Z-boson production include a factor of Bi{Z — > £~^i~) = 3.37%. 



transverse distances and is parameterized in terms of a hadronic scale Anp. For simplicity, 
we will assume that this form factor is independent of ^. The above ansatz inserts a factor 
[/hadr(2;TANp)]^ uudcr the integral over xt in fl26|) . which suppresses the region of very large 
xt values. We will employ the Gaussian model 

/S''(xtAnp) = exp (-A^p 4) (32) 

for the form factor. It was shown in j3| that the functional form of the model function only has 
a minor impact on the results, which are mainly sensitive to the value of the parameter Anp. 
Choosing Anp ~ 600 MeV shifts the position of the peak of the qx distribution for Z-boson 
production at the LHC from 3.2 GeV to 3.5 GeV and yields to a significantly better agreement 
with the data. A similar effect is seen for Tevatron data. 

In Figure O we compare the situation in Drell-Yan production of Z bosons, for which the 
characteristic scale ~ 1.75 GeV is rather low, with that in Higgs production at the LHC, 
for which g^, ^ 7. 7 GeV is safely in the perturbative domain. As expected, we find that the 
impact of hadronic effects is significantly reduced in the latter case. With Anp ~ 600 MeV, 
for instance, the peak position shifts by merely 100 MeV (from 9.1 GeV to 9.2 GeV), which 
is hardly visible on the scale of the plot. We will see in the next section that perturbative 
uncertainties are significantly larger than this effect. It is therefore safe to ignore the potential 
impact of long-distance effects for all practical purposes. 

3 Predictions for the LHC 

Having discussed the factorization of the cross section and its behavior at very small qt, we now 
present our final results for the transverse-momentum spectrum of Higgs bosons produced in 
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Figure 3: Comparison of different schemes for performing the matching to fixed-order results at 
0{as)- The bands show the resummed and matched results, the line the matching correction 
itself. 



gluon fusion at the LHC In order to obtain reliable predictions also at intermediate values, 
we match the resummed differential cross section (ITS]) to the 0{as) fixed-order result. To 
this end, we add the fixed-order cross section ctnlo to the resummed result and subtract the 
fixed-order expansion of f|T5l) so as to avoid double counting: 



da 



NNLL+NLO 



da 



NNLL 



-|- dcT] 



matching 



(io-NNLL + (c?0"NLO - da 



NNLL 



expanded to NLO 



) • (33) 



The expansion of the resummed result to 0{as) can be derived using 



2 ( 2 2 \ _ 

Q.T ^gg^ijy^ii ^2, (1ti '^hi ~2 



2T^ In ^ + 470^) 5(1 - zi) 5{l - Z2) 5,, 5gj 



(34) 



Since the hard function H(n) = C^^m^, n) ICsi—rnj^, n)f = l + 0{as), the corrections encoded 
in H{fi) do not contribute to the 0{as) expansion of the resummed cross section. Nevertheless, 
from a physical point of view one may argue that one should factor out these corrections in 
( p^ . because they are to a large extent universal. This is obvious for the prefactor C^, but it 
should also be the case for the large Sudakov logarithms and other corrections encoded in [Cgp, 



which also appear in the total cross section 2S, |25|. Factoring out the hard function i?(/x) 



before performing the matching leads to a Sudakov suppression of the matching correction at 
low Qt and an enhancement at larger transverse momentum. 

In Figure [3l we compare our results for the resummed and matched differential cross section 
obtained when the hard function is not factored out ( "naive matching" ) with the those found 
when the right-hand side of fl34p is multiplied by H{fi) ("alternative matching"). In the second 
scheme, we indeed observe the Sudakov suppression of the matching correction at low qt, but 
also a drastic reduction of the scale dependence as displayed by the width of the band. This 
is due to a cancellation of the fi dependence between the resummed result and the matching 
correction. While some reduction may be expected, such a strong effect is presumably in part 
accidental. Indeed, instead of using H{ii) with variable n we may equally well factor out the 
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Figure 4: Comparison of the resummed and matched transverse-momentum distributions of 
Higgs bosons produced at the LHC, for ^/s = 8TeV (left) and y/s = 13TeV (right). The 
default matching scheme is adopted. Note the different scales in the plots. 



hard function at the fixed default scale /j, = qx + q*- Doing so has the same qualitative effect 
on the matching correction, but as shown in the third plot the strong cancellation of scale 
dependence is not observed. To be conservative, we adopt this last choice as our "default 
matching" prescription. 

Since the matching correction for Higgs-boson production is several times larger than that 
for the Drell-Yan case (see e.g. j?!]), it would be preferable to extend the matching to the fixed- 
order cross section to two-loop order. This requires some effort, but it is possible since the 
corresponding fixed-order result is known 27H29| and has been implemented in several public 

We note that the quark beam function Iq^q{z, /i) 



codes, e.g. MCFM (30| and HNNLO [31 



Once this result is extended to the 
T5l) will be known, and 



has recently been computed to two-loop accuracy [32 
gluon channel, all two-loop ingredients for the resummed expression 
the matching should then be extended to 0{a1). 

Our final results for the resummed and matched differential cross sections for Higgs pro- 
duction at the LHC, for = 8 TeV and 13 TeV, are shown in Figure HI The shape of the 
two spectra is very similar, the main effect of the higher center-of-mass energy being an in- 
crease in the cross section by about a factor of 2. The scale uncertainty is around ±10% in 
the peak region and increases for larger g^-, as indicated in the panels below the plots. Our 
results are fully compatible with the NNLL order predictions of [4| obtained in the traditional 
CSS resummation framework [l[. The uncertainties found in that paper are slightly smaller 
in the peak region, but about a factor of 2 smaller at large qr- The reason for the reduced 
scale uncertainty is that this work implements matching to 0{a1) as well as the hard-coUinear 
two-loop corrections, which were calculated in 33 
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Figure 5: Estimate of PDF uncertainties for Higgs production at the LHC with ^/s = 8 TeV. 
The bands correspond to one standard deviation variations of the PDF sets as compiled in 
MSTW2008 NNLO (green) and NNPDF 2.3 (orange). They do not include the uncertainty 
associated with the value of a^. For NNPDF, we choose the PDF set with the MSTW default 
value as{Mz) = 0.1170. 



We have implemented the results of our calculations into a public numerical code CuTe 15 
which produces resummed and matched NNLL+NLO results not only for Higgs production, 
but also for the Drell-Yan process and the production of Z and W bosons at hadron colliders. 
The CuTe program computes cross sections with scale uncertainties and also includes different 
models for non-perturbative effects. It uses the LHAPDF interface, so it is straightforward to 
switch between different PDF sets. As an example, we show in Figure [5] the result obtained 
with NNPDF 2.3 PDF sets 3J] and compare it with our default choice. The width of the 
bands now reflects the PDF uncertainties at the level of one standard deviation. We find good 
agreement of the distributions obtained with NNPDF 2.3 and MSTW2008 NNLO, the former 
giving rise to slightly higher cross sections. 



4 Conclusions 

We have extended our previously developed formalism for the calculation of Drell-Yan cross 
sections at low and very low transverse momentum qt 0] to case of Higgs-boson production in 
gluon fusion at a hadron collider. Large Sudakov double logarithms in ttlh / Qt are resummed 
in a systematic way using RG and anomaly equations. The leading logarithms are contained 
in the hard matching coefficient Cs, which arises when the effective ggH operator is matched 
onto a two-gluon current operator in SCET. They are resummed by evolving this coefficient 
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from a hard matching scale fil ~ ~'^'h down to a characteristic scale /i ~ qt- At NLL 
order, additional large logarithms are generated by the coUinear factorization anomaly 
In transverse-position space, they give rise to an extra factor of the form (a;|,m|^)~^9s('^r'/'). 
Long-distance hadronic effects are controlled by a dynamically generated scale g*, which sets 
an upper limit on the average transverse separation in the process, (xt) ^ As long as 

the value of lies in the perturbative domain, the transverse-momentum distribution is short- 
distance dominated even for very small values of qt- Due to the difference in color factors 
{Ca vs. Cf), this scale is significantly larger in the case of Higgs production {q^: ^ 7.7 GeV) 
compared with that of Z-boson production {q^: ^ 1.75 GeV). As a result, we find that long- 
distance hadronic corrections to the shape of the transverse-momentum distribution of Higgs 
bosons produced at the LHC have a negligible effect, see Figure [21 

From a technical point of view, a novel feature of the case of Higgs production is that the 
tensor structure of the beam-jet functions B^" in ([8]) implies that the factorized cross sections 
(ITOi) and (|T2l) cannot be expressed as a simple product of two beam functions, but rather as a 
sum over products, reflecting the entanglement of the spin-0 state of two gluons [l^. However, 
we have shown that one can find a basis in which this effect arises first at NNLO in a^. 

Our numerical analysis of the transverse-momentum distribution of Higgs bosons produced 
in gluon fusion at the LHC (with ^/s = 8 and 13 TeV) presented in Section [3] gives rise to 
robust predictions, which hopefully will soon be confronted with first data from the LHC. 
We have explored different schemes for performing the matching to fixed-order perturbation 
theory and found that the various results agree within the estimated theoretical uncertainties. 
Our final predictions are presented in Figure HI with PDF uncertainties estimated in Figure O 
They suggest that at NNLL-I-NLO our predictions for the differential cross section da/dqT 
have uncertainties of order ±10% in the peak region, which seems reasonable given what is 
known about the behavior of the perturbative series for the total cross section. In order to 
reduce these uncertainties further, it would be necessary and desirable to extend our analysis 
to N^LL+NNLO, or at least to perform the matching to fixed-order theory at two-loop order. 
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